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Abstract 

Version 3 of FORM is introduced. It contains many new features that are inspired by current 
developments in the methodology of computations in quantum field theory. A number of these 
features is discussed in combination with examples. In addition the distribution contains a 
number of general purpose packages. These are described shortly. 
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1 Introduction 



It is already more than 10 years ago that the symbolic manipulation program FORM was intro- 
duced. At the time it was constructed more or less as a successor to SCHOONSCHIP H which, due 
to the limitations of its language of implementation, became harder and harder to use. In addition 
the continuous development of methods of computation in many fields of science required a rapid 
expansion of the possibilities of the symbolic programs. Admittedly programs like Reduce, Maple 
and Mathematica offer a host of possibilities combined with extensive libraries, but the speed of 
SCHOONSCHIP was always something that set it apart from the other programs. Hence the task 
was to create a program that would be at least as fast as SCHOONSCHIP and be fully portable. 
In particular it should be usable on the fastest computers that are suitable for this kind of work 
(eg. vector processors will rarely give great benefits to symbolic manipulation in general. As we 
will see this is completely different for parallel processing). This resulted in version 1 of FORM. 
Due to more modern algorithms and the benefits of a second generation approach it is actually 
somewhat faster than SCHOONSCHIP when the problems remain limited in size and much faster 
when the expressions become large. Compared to the more popular computer algebra systems the 
difference in speed and use of memory depends largely on the type of the problem. This is due 
to a completely different philosopy of design and use. Typical factors are between 10 and 100, 
but there have been rare cases for which the factor was considerably less and cases for which the 
factor was considerably larger. Such a factor translates into a timeshift with respect to when a 
given calculation can be done in practise with other programs. Over the past years the computer 
industry has given us on average every three years a factor 4 in increase in capabilities. Hence for 
FORM this translates into a shift of more than 5 years. As a consequence most of the calculations 
in quantum field theory for which the available computer resources are a consideration have been 
using FORM over the past years. 

Of course the development of FORM has continued. Version 2.0 was introduced in 1991. It 
contained a number of improvements that were inspired by the use of the program in a number of 
actual research projects. Since then regular small improvements were made, dictated by more and 
more complicated projects. After a number of years it became clear that a number of fundamental 
internal changes would enhance the capabilities very much and hence work on version 3 was started. 
Most of the code of version 2 was replaced, introducing better structured ways of dealing with 
internal objects. Even though the language of implementation is still regular C, the methods 
used are in many cases those of object oriented programming. It was of course important to do 
this in such a way that the speed and compactness of the expressions would not suffer from this. 
This has been achieved. After this revamping the new facilities were added. Some of them are 
completely new in the field of symbolic manipulation. This does not set the claim that such things 
are impossible with other programs. It just means that the FORM approach to some things is far 
more direct and often easier to program, even if some people may not like the language. 

In this paper the most important new features are described shortly. In most cases some short 
examples of code are given with the corresponding output. It should be understood that this paper 
is not a manual. Because a number of packages comes with the FORM distribution, the examples 
to illustrate the features come mainly from these packages. In addition there are examples from 
some incomplete future packages. To improve the understanding of the above, first a very short 
description of FORM is given. 

There is also a short section on a parallel version of FORM. It looks like this will be a major 
new feature that is currently still under development but results have been obtained already with 
one of the packages. The actual target is to run nearly unmodified FORM programs (the mincer 
package that has been used already was not modified at all) on a number of processors in such a 
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way that to the user it looks like a much faster sequential machine. This is very useful for -among 
others- the development cycle of big programs. 

The final section discusses the availability of FORM. In the appendix we give a short description 
of the packages with some examples of their use. 

2 A short overview of FORM 

Here a very brief overview of FORM is given. For a complete introduction one can consult the 
manual which contains a complete tutorial and a reference section. 

FORM is not an interactive system. It runs in batch. In UNIX terms one would say that 
FORM is a filter. For large programs this is the natural mode anyway. It also allows the user to 
prepare programs with a familiar editor. The editor used by the author is stedi which has a folding 
capability that is recognized by FORM. The sources of this editor are available in the FORM 
distribution as well. If everan interactive version will be constructed, it will be based on this editor. 

The language of FORM uses strong typing. This means that all objects that will be used in a 
program will have to be declared. This allows FORM to use a number of properties of such objcets 
and this adds to the speed of the program. Algebraic objects in FORM are expressions, symbols, 
functions, vectors and indices. In addition there are $-expressions, sets, numbers, preprocessor 
variables and argument field wildcards. The functions are subdivided in categories of commuting 
and noncommuting functions. Additionally there are the classes of tensors, tables and regular 
functions. 

The basic modus operandi of FORM is that variables are declared. Then one or more expressions 
are defined by giving them initial contents. After that operations are applied to one or more of 
these expressions. 

A FORM program is subdivided in modules. Each module is an individual entity of which the 
contents are compiled when it is the turn of the module to be executed. After the module has 
been executed, the compiled version is forgotten again, but depending on the type of the module, 
declarations and expressions will be kept. Hence at the level of complete modules FORM acts as an 
interpreter. The contents of expressions as they were at the start of the module will be forgotten, 
unless the user takes some special action. This is needed to prevent the disk from filling up too 
quickly. At the beginning of each module extra variables can be declared and new expressions can 
be defined. If necessary old expressions can be deleted. In principle the operations of the module 
are applied to all old and new expressions. 

The basic object for manipulation is the term. Expressions consist of terms, separated by addi- 
tion operators (or subtraction of course). This addition operator is commutative and associative. 
Function arguments can consist of subexpressions, i.e. smaller expressions which are limited in 
size. Such subexpressions also consist of terms. Terms consist of subterms. The subterms are 
separated by the multiplication operator. This multiplication operator is associative and usually 
commutative. Only in the case of noncommuting functions and still unsubstituted objects that 
could contain noncommutative functions is the multiplication order sensitive. The number zero is 
not considered to be a term. An expression that is zero contains no terms. 

FORM works its way through an expression in a sequential manner. It starts with the first term, 
applies all statements of a module to it and puts the results away for sorting (normal ordering of 
the complete expression is called sorting here. It does involve adding coefficients of terms that are 
identical with the exception of their prefactor) . Partial sorting may occur when buffers become too 
full. After all terms of an expression have been treated this way the sorting is completed. Once all 
expressions have been treated this way the new versions of the expressions replace the old versions. 
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The actions inside a module can be either predefined operations or substitutions that may in- 
volve rather complicated patterns involving wildcarding. Wildcarding is the use of generic variables 
that represent a whole category of objects, both in the left hand side (the pattern) and in the right 
hand side (the replacement expression). 

A number of special functions is available to make the use of the program easier and often much 
faster. These involve amoung others a sum function, a number of combinatoric functions and a 
number of numerical functions. In addition there are special tensors like the Kronecker delta and 
the Levi Civita tensor. And there is much more. 

For the use in particle physics there are Dirac gamma matrices and trace operations. It should 
be noted however that this is nowadays a rather unimportant part of the program as there exist 
many facilities to emulate these objects and operations (and similar operations that can be much 
more complicated). 

Results of a calculation can be presented in many different ways. One way is to print them to 
regular output. They can be written to file in internal notation for use in future FORM programs. 
Or they can be written to special files, etc. 

There is a powerful preprocessor that can manipulate the input at the character level. It is fully 
adapted to the requirements of symbolic manipulation. There are now also some mechanisms for 
the feedback of algebraic results to the preprocessor. This facilitates the possibility to have FORM 
programs write parts of themselves. 

A number of buffers in FORM have to be allocated at startup. The user can influence these 
allocations with a special setup file. Ordinarily this is considered an advanced feature. 

3 New features 

For the description of the new features it is assumed that the user has some experience with the 
older versions. The distribution contains a good tutorial. The interested reader without experience 
could look at the many examples in it. Here we concentrate on the new concepts. For a complete 
description and the exact syntax and for the many additional new features we refer to the new 
manual. 

3.1 Variable administration 

The complete internal adminstration of variables has been made dynamic. This means that the 
old restrictions on the number of variables have been lifted. The setup file is no longer needed 
when more than 100 variables of the same type are asked for. The limit lies now at 6000 for 32-bit 
systems and 10 8 for 64-bit systems. 

Similarly the compiler buffer, the input buffer and the various preprocessor buffers are arranged 
dynamically. The reason that in the older versions these things were set up in a static way was due 
to the size of the computer memories in the late 80's. 

3.2 Preprocessor variables 

Preprocessor variables can be used 'recursively'. The old notation to enclose a preprocessor variable 
in quotes has been changed into a notation in which the variable is enclosed in a backquote-quote 
pair as if these are brackets. This way they can be nested as in 

#do itabs = 1,'size' 
#ifndef 'STABLE' itabs 'HFILE' 
#call table 'itabs' 
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#endif 
#enddo 



In the summer package there are 8 files table'itabs'.prc all having something like: 

#procedure table5 
#ifndef ' STABLE5HFILE' 
#define STABLE5HFILE "1" 

lots of tables 
#endif 

#endprocedure 

This way the procedures are not called twice. Of course the contents would not be used twice 
anyway, but some of them are rather lengthy and FORM would have to read the whole procedure 
to find the #endif. That is avoided this way. 

To make the preprocessor variables even more powerful there are also postincrement and post- 
decrement operators for preprocessor variables as in the C language: 

#define k "1" 
Local F'k++> = . . . ; 
Local F'k++> = . . . ; 
Local F'k++' = . . . ; 
Local F'k++' = . . . ; 

would result in the definitions of Fl to F4. In the end k would have the value 5. This is particularly 
handy when the program is generating expressions and we do not know in advance how many. 

3.3 Preprocessor abbreviations 

The notation has been cleaned up a bit by introducing a new operator for the compilation stages. 
This operator is indicated by three dots as in 

Symbols xl , . . . , xlOO ; 
Local FaclO = 1* . . . *10; 

Actually there is an even more general notation involving <>: 

id f (<pl?,ml?,il?>, . . . ,<p6?,m6?,i6?>) = 

Here the <> enclose a pattern. FORM will look for the differences in the begin pattern and the 
end pattern, calculate their numerical differences and if they all have the same absolute value we 
can make a 'running pattern' and make the substitution. Hence 

id f (<pl,m4>, . . . ,<p4,ml>) = 

would evaluate into 

id f (pl,m4,p2,m3,p3,m2,p4,ml) = 

This facility cannot be used recursively but already it can clean up the code considerably when 
compared to how one would have to do this with do loops in the old versions of FORM. It becomes 
really powerful in combination with preprocessor variables: 
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#do j = 1, 'MAX'-l 

id e_(il?,...,i'j'?) = 

sum_(k, 1 , 'MAX' , e_ (il , . . . , i ' j ' ,k) *f ( ' j '+1 ,k) ) ; 

id f('j\k?) = tabCj'.k); 

. sort : step ' j ' ; 
#enddo 

id e_(l, . . . , 'MAX') = 1; 

This would be the central step in computing the determinant of a matrix that has been defined in 
the two dimensional table named tab. We use an intermediate function f to make sure that first 
the Levi-Civita tensors with two indentical indices are killed before table substitution is attempted. 
This is a very efficient algorithm for sparse matrices, especially when the matrices have been ordered 
in such a way that the first columns have the largest number of zeroes. 

3.4 $ variables 

These variables are at the same time small algebraic expressions and character strings. They can 
be used as either. When they are used and they are enclosed between a backquote-quote pair, the 
compiler will make their character representation and substitute this like a preprocessor variable. 
This use as a preprocessor variable even includes the possibility of using the postincrement or 
postdecrement operators. If they are not enclosed between the backquote-quote pair the contents 
are seen as an algebraic expression to be used during execution time. This can lead to very flexible 
code. If they are to be given an initial value, this can be done with the character # before their 
name to have the preprocessor provide the value (hence this is done once during compilation time) , 
or without the # to indicate that this value should be provided during execution, each time a term 
passes this position in the code. There are additional methods to provide the $-variables with a 
value that is the match of a wildcard variable. Let us have a look at an example: 

#$max = -100; 

if ( count(x,l) > $max ) $max = count_(x,l); 
. sort 

#do i = 1 , ' $max ' 
#enddo 

Here we give $max the initial value -100. During execution we look for the maximum power of 
the variable x. Once the module has finished this maximum should be indeed the value of $max, 
provided it is at least -100. At this point we can use this value as a parameter in a do loop. Note 
that the .sort is essential, because if we were to omit it, the do loop would have the bounds 1,-100 
because the do loop instruction would be translated before any algebraic manipulation would have 
taken place (first the complete module is translated, then execution takes place). 

It is also possible to manipulate the contents of the $-variables, provided these do make sense. 

The contents of these $- variables are kept in memory. This means that the user should not try 
to see them as a substitute for regular expressions and have some of them grow to a large size. 
It would work, but performance might suffer when the physical memory is exhausted in the same 
way that Mathematica and Maple become rather disagreeable when they need more memory than 
is available. 
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3.5 Hide 



When there are several expressions active at the same time, sometimes there is the need to tem- 
porarily put some of them 'out of the way'. FORM has several mechanisms for this. The Skip 
statement would just skip the indicated expression(s) by copying them directly from the input to 
the output of the module. This is not always efficient. In addition one has to specify this in each 
module seperately. A more basic solution is to copy the expression(s) to a special location where 
they remain suspended until called back. This is done with the Hide and Unhide statements. The 
hidden expressions can still be used in the right hand side of expressions in the same way as fully 
active expressions can be used. This is a faster mechanism than the one that is used with stored 
expressions. 

. sort 

Hide F1,F2; 
Local F4 = Fl+f3; 
. sort 

more modules 

. sort 

Unhide F1,F2; 
id x = y+1; 

Here Fl and F2 are not operated upon for a few modules. Then they are called back and they are 
immediately active again. This means that the substitution of x also holds for the terms of Fl and 
F2. 

There is another mechanism when one or a few statements should be applied only to a limited 
number of expressions. This concerns a new option of the if statement: 

if ( expression(G) == ) ; 

statements 
endif ; 

In this case the statements are applied to the terms of all active expressions except for the terms 
of the expression G. 



3.6 #write 

This is a preprocessor controled instruction that can write either to regular output, to the log file 
or to any named file. It has a format string as in the printf function of the C language. It specifies 
what should be printed. There are some examples in which complete procedure files are constructed 
this way for later use in the same program. 



3.7 Print "" 



The new feature of the print statement is a mode like printf in the C language that can print 
information during execution. It can handle the printing of messages and even individual terms 
or $-variables. This gives it great value as a debugging tool. A good example is given with the 
Split Arg statement in 3.11. 
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3.8 Table manipulations 

This facility involves a statement, a special function and a special print statement. The table, 
function can be used to load the contents of a table into an expression. The table indices are 
represented by powers of a given set of symbols. The fillexpression statement does the opposite. It 
takes an expression that is bracketed in a number of symbols and loads the contents into a table. 
The powers of the symbols indicate the table indices and the contents of the brackets form the 
corresponding table elements. The printtable statement writes a table to file in the form of fill 
statements. This file can be used in a future program to fill the table. Together these facilities 
allow the dynamic extension of tables in a program. If a table element that is needed has not yet 
been defined, one could compute it on the spot and add it to the table. The results can be written 
in such a way that in future programs the same element does not have to be recalculated. 

3.9 Symmetry properties 

Functions can have properties of symmetry. These properties hold for the set of all arguments. The 
functions or tensors can be symmetric, antisymmetric, cyclesymmetric and r ever secy clesymmetric. 
This last means that the function is invariant both under cyclic permutations and a complete 
reversal of all arguments. These properties are put in the declaration of the function or tensor. Each 
time some arguments of the function are changed, FORM will check whether permutations have to 
be made. There are some restrictions. FORM will not accept argument field wildcards (wildcards 
that match zero, one or more function arguments) in patterns in which such a wildcard occurs 
inside a symmetric or antisymmetric function. The pattern matching would become extremely 
slow when no match occurs (all possibilities have to be searched and that involves n! possibilities 
for n arguments). There are tricks to circumvent this in most cases. 

3.10 Topology 

In combination with the symmetry properties FORM has now a command to look for loop structures 
in argument contractions. What this means is that if we have a network in which v is a vertex 
function and its arguments are indices and connections are represented by two vertices having 
identical indices, FORM can find the smallest loop in such a structure and make replacements 
based on this. There are some examples in the Color package in which the vertex function is the 
antisymmetric tensor f with three indices. In the program 

#- 

Indices il , . . . , i9; 

CF f (antisymmetric) ,ff ; 

Local F = f (il,i2,i3)*f (i2,i4,i5)*f (i3,i5,i6)* 
f (i4,i7,i8)*f (i6,i7,i9)*f (il,i8,i9) ; 
ReplaceLoop , f , arguments=3 , loopsize=all , outf un=f f ; 
Print +f; 
. end 

the replaceloop command takes the smallest loop, removes the loop indices and puts the remainder 
in the output function ff which should be at least cyclesymmetric (at least the action we take here 
implies it). If there is more than one possibility FORM takes the first one it encounters. The result 
of this program is 

F = f (il,i8,i9)*f (i4,i7,i8)*f (i6,i7,i9)*ff (il,i6,i4) ; 
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Only a single loop has been replaced. If one would like to replace all loops one should put a 
repeat/endrepeat combination around the replaceloop statement. This would result in 

F = - ff (il,i4,i6)*ff (il,i6,i4) ; 

If such a command would not exist, the color package would be much more restricted because it 
would need a statement for each loop (loops with one vertex, two vertices, three vertices etc.) Also 
the corresponding pattern matching would be much slower. In the crucial phases of the program 
this made a factor 10 difference in speed and it also made the code much shorter. 

3.11 SplitArg 

This is a rather useful tool for investigating arguments of functions. In its basic form it splits a 
multiterm argument into a number of arguments, each with a single term. There are variations 
for faster operation when only a single term should be taken from the multiterm argument. Take 
for instance the following example in which we want to select all functions den that contain the 
symbol j 1 in its argument. 

SplitArg, ( (j 1) ), den; * take all terms in jl separately 

id den(jl) = denl(0,jl); * a single term does not get split 

* we may need more of these statements 

id den(x?,xl?) = denl (x/xl*j 1 , j 1) /xl*j 1 ; * normalize the jl argument. 
* 

* all den with two arguments have now just jl without factor. 
* 

repeat id denl(x?, jl)*jl = l-x*denl(x, jl) ; 
repeat id denl (xl? ! {x2?} , j 1) *denl (x2? ! {xl?} , j 1) = 

(denl (xl , j 1) -denl (x2 , j 1) ) *den(x2-xl) ; 
id den(x?number_) = 1/x; 

Code of this nature occurs frequently in the summer package. It splits the fractions, being care- 
ful not to divide by zero (in some cases one has to be far more careful and involve theta_ and 
delta, functions). Let us see what this code does. We use the following program (this is also a 
demonstration of the new print statement capabilities): 

#- 

S jl,x,xl,x2; 
CF den, denl; 
Off Statistics; 

L F = den(jl)*den(2+jl)*den(3-2*jl) ; 

Print +f "<1> %t"; 
SplitArg, ((jl)) , den; 

Print +f "<2> %t"; 
id den(jl) = denl(0,jl); 

Print +f "<3> 7„t"; 
id den(x?,xl?) = denl (x/xl* j 1 , j 1) /xl*j 1 ; 

Print +f "<4> %t"; 
repeat id denl(x?, jl)*jl = l-x*denl(x, jl) ; 

Print +f "<5> %t"; 
repeat ; 
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id denl(xl?!{x2?}, jl)*denl(x2?!{xl?},jl) = 
(denl (xl , j 1) -denl (x2 , j 1) ) *den(x2-xl) ; 
Print +f "<6> %t"; 
endrepeat ; 

id den(x?number_) = 1/x; 

Print +f "<7> %t"; 
Print +f; 
. end 

which results in 
#- 

<1> + den(jl)*den(2 + jl)*den(3 - 2* j 1) 

<2> + den(2, jl)*den(3, - 2*j 1) *den( j 1) 

<3> + den(2, jl)*den(3, - 2*j 1) *denl (0 , j 1) 

<4> - l/2*denl(0, jl)*denl(2, jl)*denl( - 3/2, jl) 

<5> - l/2*denl(0, jl)*denl(2, jl)*denl( - 3/2, jl) 

<6> - l/2*den(2)*denl(0, jl)*denl( - 3/2, jl) 

<6> - l/2*den(2)*den( - 3/2)*denl(0, jl) 

<6> - l/2*den(2)*den( - 3/2)*denl(0, jl) 

<7> + l/6*denl(0, jl) 

<6> + l/2*den(2)*den( - 3/2)*denl( - 3/2, jl) 
<6> + l/2*den(2)*den( - 3/2)*denl( - 3/2, jl) 
<7> - l/6*denl( - 3/2, jl) 

<6> + l/2*den(2)*denl(2, jl)*denl( - 3/2, jl) 
<6> + l/2*den(2)*den( - 7/2)*denl(2, jl) 
<6> + l/2*den(2)*den( - 7/2)*denl(2, jl) 
<7> - l/14*denl(2, jl) 

<6> - l/2*den(2)*den( - 7/2)*denl( - 3/2, jl) 
<6> - l/2*den(2)*den( - 7/2)*denl( - 3/2, jl) 
<7> + l/14*denl( - 3/2, jl) 

F = 

- 2/21*denl( - 3/2, jl) + l/6*denl (0 , j 1) - l/14*denl (2 , j 1) ; 

Note that the termwise option of the print statement gives us a trace of what happens and in what 
order. 

3.12 The term environment 

Sometimes one needs to work things out for just a few terms, but an intermediate sort for the 
results of these manipulations is needed in order to have things cancel. Or one would like to do 
even fancier things like trying whether the contents of certain brackets can be factorized by a given 
expression. If so, one would like to make that factorization. The term environment defines the 
current term temporarily as a small expression. Because this is done during the processing of the 
terms, sorting cannot be done with the .sort instruction which works for the whole expression. This 
means that we need a separate sort statement for all these small expressions. It should however be 
noted that the complete term environment consisting of the term statement, the endterm statement, 
and the intermediate regular statements and sort statements are all part of the same module and 
hence are all compiled before execution starts. Basically we use the same mechanism of operation 
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in the argument environment. There function arguments can be operated upon and when the 
endargument statement is encountered the results are sorted and put back inside the function. 
The sort statement is hence equivalent to an endargument statement followed by a new argument 
statement. This last construction has the limitation that the result has 'to fit' and hence is limited 
to MaxTermSize or less. The Term environment does not suffer from this restriction because we 
do not have to keep the results inside a function argument and also in the intermediate stages even 
the disk can be used and the expression may be quite large. 

#- 

S x,y, [x+1] , [x+2] ,xl,x2; 
CF acc; 

L F = y*(x+l)~2*(x+2)+y~2*(x+l)*x+y~3*(x+2)*x+y~4*x~2; 
B y; 
Print +f; 
. sort 

Collect, acc; 
Term; 

$minl = 1000; 

$min2 = 1000; 

id acc(x?) = x; 

id x = xl-1; 

sort ; 

if ( count (xl,l) < $minl ) $minl = count_ (xl , 1) ; 
sort ; 

Multiply ([x+l]/xl)~$minl; 
id xl = x+1; 
id x = x2-2; 
sort ; 

if ( count (x2,l) < $min2 ) $min2 = count. (x2, 1) ; 
sort ; 

Multiply ( [x+2] /x2) ~$min2 ; 

id x2 = x+2; 
EndTerm ; 
Print +f; 
. end 

which gives 
#- 

F = 

+ y * ( 2 + 5*x + 4*x~2 + x~3 ) 

+ y~2 * ( x + x~2 ) 

+ y~3 * ( 2*x + x~2 ) 

+ y~4 * ( x"2 ) ; 

F = 

x*y~2*[x+l] + x*y~3*[x+2] + x~2*y~4 + y* [x+1] "2* [x+2] ; 
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Of course this works only if one has a suspicion which factors are involved and things would go 
much better if FORM would be able to deal with factorizations by itself. This is considered for the 
future. For now this provides at least a partial solution. 

3.13 Bracket indexing 

In the old versions FORM would look for the contents of brackets inside an expression in a sequential 
manner. For relatively small expressions in computers with not much memory this was still not 
optimal, but also not a complete disaster. With the growing of the computers and hence the 
expressions the inefficiency of this method became a major factor as it is an algorithm that has 
a linear dependence of the size of the expression. Hence now the brackets can be indexed if the 
user asks for it. The index will of course take space of its own, and this space can approach (in 
principle) the size of the expression if the brackets have very little contents. This could result in 
bad problems again. Therefore there is a maximum to the size of the index (can be changed in the 
setup file) and when this maximum is reached FORM will start skipping brackets in the index. If 
such a skipped bracket is asked for, FORM will do a linear search, starting from the nearest bracket 
in front of it. The result is still a speedup by a big factor. The brackets in the index are located 
by a tree-like binary search which gives a logarithmic dependence of the number of brackets in the 
index. The degrading factor when brackets are skipped is proportional to the size of the expression 
divided by the number of brackets in the index. Hence in that case the improvement is proportional 
to a factor B/ln(B) with B the number of brackets in the index. Usually however there should 
be no problem in fitting the entire index inside the memory. Because this use of brackets is not 
very common and the making of the index does cost resources, the default is to not make an index. 
Here we have a simple example 

#- 

Symbols xl , . . . ,xlO; 

Local F = (xl+. . .+xlO) "10; 

Bracket+ xl; 

. sort 

Drop F; 

#do i = 0,10 

Local F'i' = F[xl~'i>] ; 

#enddo 

. end 

When we run this program the second module takes 1.59 sec on a pentium 300. Without the + in 
Bracket-)- this step takes 2.84 seconds and this is not a very big expression yet. The program that 
solves the equations for the values in N — > oo in the summer package had an improvement factor 
of more than 100. 

3.14 Error messages 

The error messages have been changed completely. The most noticeable change is the indication of 
where the error occurs when the input listing is off. This mode is standard for procedures. There are 
additional warnings and there is a way to let FORM print its interpretation of the input. Together 
with the already existing printing of the name lists this provides extra tools for the debugging of 
FORM programs. In this example we have made two errors. It is the program of the bracket index. 
The name of the file is exl.frm. Because we had turned off the listing, the line number information 
is rather crucial: 
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#- 

Symbols xl,...,xlO; 

Local F = (xl+. . .+xl0)~~10; 

Bracket+ xl; 

. sort 

Drop F; 

#do i = 0,10 

Local F'i' = FF[xl~ f i'] ; 

#enddo 

. end 

When this program is run we obtain: 
#- 

exl.frm Line 3 — > Illegal position for operator: "10 
exl.frm Line 8 — > Undeclared variable FF 
++++Errors in Loop 

Because nearly all error messages have been changed and actually the complete compiler has 
been replaced, it is possible that some error messages either do not make complete sense or are too 
vague. The author is open for suggestions. 

4 The parallel version 

Because of the way FORM processes expressions it should be possible to achieve a high degree of 
parallelization for the complete program. After an initial study at the special computer at FNAL in 
1991 this assumption was verified. It took however till the late 90's before a more serious attempt 
was made. This was made possible in the framework of a large phenomenology project centered 
at the university of Karlsruhe (DFG Contract Ku 502/8-1). Two people worked part time on this 
(Albert Retey and Denny Fliegner) and after struggling for a long time with the hardware and the 
parallelization library software they did indeed manage to obtain a working version of FORM that 
can run the Mincer package on an eight processor machine with a speed that ranges from 3 to 5 
times the speed of a single processor. In addition many tests were made with other architectures 0]. 

It should be noted that thus far the parallelization has been rather crude. The terms of the 
input expression in a module are distributed over the various processors. These do their work, and 
they sort their results. After this the outputs are collected by the master processor for a final merge 
and the output is written to disk. It is possible to be far more refined about this and make a much 
better load balancing by interfering, if needed, in the main expansion tree. The main problem 
however is the speed of transfer of operation. Thus far use has been made of the MPI protocol 
(and in tests of much more of course). In the future it is expected that the use of threads will 
increase the data transfer rate considerably. It is of course this transfer rate that puts a limit to the 
eventual speed that can be obtained. A given transfer speed in combination with slow processors 
will lead to saturation for a much larger number of processors than in the case that fast processors 
are used. In a shared memory machine with threads however the transfer rate would be mostly 
a function of the processor speed and hence much is expected of this model. Unfortunately the 
implementation is not yet complete. 

At the moment the only readily available version is the one for a Compaq Alpha server (with 
8 Alpha 21264 processors). It is expected that more will follow. A time prognosis cannot be given 
yet. 
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Of course it is not always possible to run code in parallel. Especially the new $-variables could 
cause problems because their eventual contents depend on the order in which terms passed. If each 
processor keeps its own version of a $-variable the processors may not have identical values at the 
end of the module. Take the simple code 

#$max = -100; 

if ( count (x,l) > $max ) $max = count_(x,l); 

which tries to find the maximum power of x. The various processors may have a different value in the 
end. Hence FORM should not parallelize modules in which $- variables obtain values during running 
time, unless the user can provide some extra information. This can be done in the moduleoption 
statement. In that statement one can tell FORM that the final value of the $- variable is either a 
sum, a maximum, a minimum or completely unimportant. In such cases FORM knows what to do 
with the results and the module can be run in parallel after all. It should be realized that this is 
all very new. Hence it could well be that many new facilities will be developed here and that there 
may be more ways in which future versions of FORM can decide about parallelization. 

Any statement that concerns the parallelization will be completely ignored by a sequential 
version of FORM. This means that when a program is being prepared one could or should include 
already some of the indications for the parallel version. This way the upgrading to a parallel 
machine would be rather smooth. 

5 Availability 

Whereas version 2 of FORM has been commercial, it was judged that the scientific community as 
well as the reputation of the author would benefit more from a free distribution. Hence, starting 
version 3, FORM will be freely available again. For the moment this will involve a number of 
binary executables of the program. Amoung these will be at least executables for LINUX on PC 
architectures and for alpha processors running UNIX. More binaries will become available in the 
near future. The manual will be available both as a postscript file and as a .pdf file. This will 
allow the average user to obtain the system and print out its manual as well as view it on the 
computerscreen. 

Of course it is hoped that this development will not only improve the visibility and popularity 
of FORM, but will also lead to more possibilities to improve FORM itself. The potential for this 
depends of course on people in charge of jobs knowing about the importance of FORM in the 
science community. Hence it is expected that people who use FORM for scientific publications 
will refer to the current paper. This way the use of FORM can be 'measured'. The author hopes 
that the free distribution in combination with a proven popularity will lead to more people being 
involved in the development of FORM itself. There is still a whole list of potential improvements 
waiting for implementation. Undoubtedly many users will have their own suggestions. Also more 
packages would be appreciated. 

The FORM distribution contains: 

• Executables for a variety of computers. 

• The conv2to3 program with its sources. This program converts old FORM code into new 
FORM code. 

• A postscript and a pdf version of the manual and a tutorial by Andre Heck. 

• A number of packages. Currently these are 'summer', 'harmpok, 'meltran', 'color' and 'min- 
cer'. 
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• The sources of the stedi editor. This editor works on BSD 4.0 systems and systems that still 
have an old compatibility mode. It also works under LINUX. There is a manual with it. 
Originally it was an editor for the Atari ST computers, made to work as well under MS-DOS. 
Then it was ported to BSD UNIX where it did run inside terminal windows. Under LINUX 
and X-windows it has recovered most of its old ease of use. 

• The sources of the minos database program for organizing large calculations. There is a small 
manual with it. The program should work on nearly any system. The database files should 
be system independent. 

• The axodraw system for including simple drawings and figures inside DTpXfiles. This includes 
the manual. 



The files can be obtained from the www pages in http:/ / www.nikhef.nl/~form . 
Disclaimer: The files are provided without guarantee whatsoever. If some commands do not work 
as mentioned in the manual or as the user thinks they should work, the user should submit a bug 
report. This may or may not be reacted to depending on circumstances. 

Acknowledgements: The author is thankful for many comments during the development of 
FORM. Many people contributed in the form of bug reports and useful suggestions. Some people 
however stand out. Denny Fliegner and Albert Retey for developing the parallel version of FORM 
(with support of DFG Contract Ku 502/8-1 (DFG-Forschergruppe 'Quantenfeldtheorie, Comput- 
eralgebra und Monte Carlo Simulationen') , Geert Jan van Oldenborgh for being the main guinea 
pig during the early stages of the project, Andre Heck for writing a good quality tutorial, Walter 
Hoogland for allowing me to spend so much time on the project in its early stages when success 
had not been proven yet. 



A The packages 

A few packages are provided with the FORM distribution. These packages are useful not only 
for their field of applicability. They also serve as examples of FORM packages and programming 
techniques. In this paper we took most of the examples from them. Here we give a very short 
description of what these packages are good for. More information is in a number of papers and in 
the future there may be some manuals. 



A.l Summer 

This packages deals with harmonic sums. A full description of these sums is given in ref. . We 
give a short overview. The harmonic series is defined by 

n 1 

Sm(n) = £- (!) 
i=i 1 

n f-lV 

S-m{n) = E^T- ( 2 ) 
i=l 1 

in which m > 0. One can define higher harmonic series by 

n 1 

s m,n,-,j P ( n ) = X TT^i. • •,„.(') ( 3 ) 
1=1 L 

n 1-1 V 

S-m, h ,..,i P {n) = J2 K -^S hj .., jp (i) (4) 

1=1 
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with the same conditions on m. The m and the ji are referred to as the indices of the harmonic 
series. In the program we will use Sj lr ..j p (n) — > S(R(jl, . . . , jp) ,n). 

There is an alternative notation in which the indices are one of the values 0,1,-1. A zero indicates 
that actually one should be added to the absolute value of the nonzero index to the right of the 
zero as in: 

50,1,0,0,-1,1 ( n ) = S 2,-3,iW (5) 

The weight of a sum is sum of the absolute value of its indices. In the notation with the 0,1,-1 
it is equal to the number of indices. The harmonic sums with the same argument form a (weight 
preserving) algebra: The product of two sums with weights w\ and u>2 respectively is a sum of 
terms, each with a single sum of weight W\ + W2- The routine 'basis' in the package implements 
this property. Example: 

#- 

#include summer6.h 
.global 

Local F = S(R(2,3) ,N)*S(R(-1,2) ,N) ; 
#call basis (S) 
Print +f; 
. end 

This program gives the output 

F = - S(R(-3,2,3),N) - S(R(-3,3,2) ,N) + S(R(-3,5),N) 

+ 2*S(R(-1,2,2,3) ,N) + S(R(-1,2,3,2) ,N) - S(R(-1,2,5) ,N) 

- S(R(-1,4,3) ,N) - S(R(2,-4,2) ,N) + S(R(2,-1,2,3) ,N) 

+ S(R(2,-1,3,2) ,N) - S(R(2,-1,5) ,N) + S(R(2,3,-1,2) ,N) ; 

This defines a standard basis. A large number of sums involving harmonic sums and denominators 
can be done to any level of complexity. We give here a relatively simple example: 

#- 

#include summer6.h 
.global 

Local F = sum(j ,1,N)*S(R(2,3) ,0, j)*S(R(-l,2) ,N,-j)*den(0,j) ; 
#call summit () 
Print +f; 
. end 

which gives the output 

Time = 0.15 sec Generated terms = 104 

F Terms in output = 78 

Bytes used = 2938 

F = - S(R(-l,-3,2,3) ,N) + S(R(-l,-2,-2,-l,3) ,N) + S(R(-l,-2,-l,-2,3) ,N) 

- S(R(-l,-2,2,-4) ,N) + S(R(-1, -2,2, -3,-1) ,N) + S(R(-l,-2,2,-2,-2) ,N) 

+ S(R(l,2,-l,-3,-2) ,N) + S(R(l,2,-l,2,-3) ,N) + 2*S(R(1 ,2,-1 ,3,-2) ,N) 
+ 3*S(R(1, 2, -1,4,-1) ,N) - 4*S(R(1,2,-1,5) ,N) + S(R(1,2,3,-1,2) ,N) ; 
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When iV — > oo there are relations in addition to the regular algebraic relations. They allow all 
the sums in infinity to be expressed in terms of a rather small number of trancendental numbers 
of the Euler Zagier type 0. In order to obtain these expressions one has to solve for each weight 
a number of equations. This number grows exponentially with the weight. The system has been 
solved exactly up to weight 9 and the file with the solutions is part of the distribution. It is about 
20 Mbytes. In the file 'summer6.h' only the solutions up to weight 6 have been included. The 
programs that solve these systems of equations are also part of the distribution. In the case of 
weight 9 there were more than 40000 equations in 13122 objects. The divergent sums were part of 
this because they can all be expressed properly into combinations of the basic divergence S\(oo). 
This database of sums is rather important for whole categories of integrals as we will see in the 
harmpol and meltran packages. It should be noted that at the moment of writing, no constructive 
algerithms are known to determine these constants one by one. The only 'algorithm' that comes 
close is the numerical evualuation of each individual sum, making an ansatz about which constants 
form a minimal set and then fitting all sums with a special numerical program [Q]. This has been 
done for w = 10 and in a restricted way for some higher values of w [||. 

A. 2 Harmpol 

This is a package for Harmonic Polylogarithms as defined in ref. || . A short description is in order 
here. 

The weight one harmonic polylogarithms are given by: 

H(0;x) = Ins, 

px A T I 

H(l;x) = / - = -ln(l-s), 

' 1 — x' 

dx' 

1 + x' 







dx' 

H(-l;s) = / — — = ln(l + x). (6) 



For the higher weights we first define 



f(0;x) 



1 

X 



1 



Then 



1 — x 

/(-i;*) = (7) 



H{5 w ;x) = —,W°x, (8) 



with ftu, an array with w zeroes while if the sequence fh w of w elements 0,1,-1 is not equal to 6 W 

rx 

R(rh w ;x) = / dx' f(m w \x) H(m UI _i;x / ) . (9) 
J o 

The weight of the harmonic polylogarithms is the number of indices. The argument x can in 
principle take any complex value, but depending on the indices there may be a cut for real values 
> 1 or < 0. These cuts are due to powers of ln(l — x) or ln(x). Such powers can be extracted, 
resulting in combinations of these logarithms and H-functions without these cuts. This will be 
needed for the Mellin transforms. It should be noted that we can use also the notation of the 
harmonic sums for the indices (with values beyond the set 1,0,-1). 
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Also the harmonic polylogarithms form a weight preserving algebra. In addition there are extra 
relations for x = 1. One has to be rather careful with the divergent elements at x = 1. The 
powerseries expansion of H contains harmonic sums. The algebraic relations for general x are 
related to the N — > oo relations of the harmonic sums and the extra relations in x = 1 are related 
to the general algebraic relations for the harmonic sums. 

The transformations x — > 1 — x, x — > 1/x and x — > give some interesting relations that 
allow for a reasonably fast numerical evaluation of the function over most of the complex plane 
(excepting the regions around ±i). 

The values in x = 1 can be expressed in terms of the harmonic sums in N — > oo. Because a 
weight w H-function is an integral over a denominator and a weight w — 1 H-function this gives 
us whole categories of integrals that were till now extremely hard to solve. The results will be in 
terms of a limited number of trancendental numbers, but the relations for x — > 1 — x at x = \ will 
allow them to be expressed in terms of relatively fast converging series. Hence these numbers can 
be evaluated to any precision when the need arises. 

Because the harmpol package uses a number of the routines of the summer package it will 
automatically load this package. Hence once harmpol.h is included, so is summer6.h. Example: 

#- 

#include harmpol . h 
Off statistics; 
. global 

Local F = H(R(1,0,1) ,x)*H(R(-l , 1 ,-1) ,x) ; 
#call hbasis(H,x) 

repeat id H(R(?a,n? ! {1 ,0, -1} , ?b) ,x?) = H(R(?a,0,n-sig_(n) ,?b) ,x) ; 
. sort 

On statistics; 
Print +f; 
. end 

gives: 

#- 

Time = 0.78 sec Generated terms = 14 

F Terms in output = 14 

Bytes used = 468 

F = 

H(R(-l,l,-l,l,0,l),x) + H(R(-l,l,0,l,-l,l),x) + 2*H(R(-l,l,0,l,l,-l),x) 

+ 2*H(R(-1, 1,1,-1, 0,1) ,x) + 2*H(R(-1, 1,1, 0,-1,1) ,x) + 2*H(R(-1 , 1 , 1 , , 1 , 
-l),x) + H(R(1, -1,0, 1,-1,1) ,x) + 2*H(R(1, -1,0, 1,1,-1) ,x) + H(R(1,-1, 1,-1 
,0,1), x) + H(R(l,-l,l,0,-l,l),x) + H(R(l,-l,l,0,l,-l),x) + H(R(1 ,0,-1 , 1 , 
-1,1), x) + 2*H(R(l,0,-l,l,l,-l),x) + H(R(l,0,l,-l,l,-l),x); 

The conversion of combinations of denominators, logarithms and polylogarithms in terms of H- 
functions is usually a straightforward exercise in integration by parts, provided that it can be done 
at all. We give one example of an integral 

#- 

#include harmpol.h 
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#call htables(6) 

CF int ; 

Off Statistics; 

.global 

Local F = ln_(l+x)~2*ln_(l-x)~2*ln_(x)/x*int(x,0,l) ; 
id ln_(l+x) = H(R(-l),x); 
id ln_(l-x) = -H(R(l),x); 
id ln_(x) = H(R(0),x); 
#call hbasis(H,x) 

id H(R(?a) ,x)/x*int(x,0, 1) = H(R(0,?a) , l)-H(R(0,?a) ,0) ; 

id H(R(?a,n?{l,-l},?b) ,0) =0; * Provided there is a nonzero index 

#do i = 1,6 

id H(R(nl?, . . . ,n'i'?) ,1) = Htab ' i ' (nl , . . . ,n< i ' ) ; 

#enddo 

. sort 

On Statistics; 
Print +f; 
. end 

which gives 
#- 



Time = 0.61 sec Generated terms = 10 

F Terms in output = 10 

Bytes used = 172 

F = - l/2*z2*ln2~4 - 129/140*z2~3 + 7/6*z3*ln2~3 - 37/16*z3~2 
- 31/8*z5*ln2 + 8*ln2*li5half + 4*ln2~2*li4half 
+ l/9*ln2~6 + 8*li6half + 2*s6; 

The constants are explained both in the paper and the htable files. These integrals are not easy by 
conventional techniques. 



A. 3 Meltran 

This package does the Mellin transforms and inverse Mellin transforms between the harmonic 
polylogarithms and the harmonic sums. It allows calculations to be done in either space and then 
converted to the other, whatever is easiest. The definition of a Mellin transform is: 



M{f{x),N) = J 1 dxx N f(x) 



^(l-X) + , JO l-X 

M ^r\ N) = (*"'<" (10) 

(1— x) + Jo 1-x 

in which the function / is supposed to be finite for x = 1 when the factor 1/(1 — x)+ is present. 
Hence we see that the first step in the Mellin transformation of a H-function is the extraction of 
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potential powers of ln(l — x). The Mellin transform of a H-function, divided by either 1 — x or 1 + x 
can be solved recursively and the result is a combination of harmonic sums in the Mellin parameter 
N. There is a one to one correspondence between the Mellin transform of H-functions of weight w, 
divided by either 1 — x or 1 + x and harmonic sums of weight w + 1 (the Mellin transform can also 
contain sums of a lower weight, but of the highest weight there is only a single term). Hence the 
inverse Mellin transform can also be obtained rather easily. The only complicating factor is that 
there will be terms in x = 1 or N — > oo. Currently these can only be combined into a 'minimal' 
representation if the weights are at most 9. 

The mellin package is part of the harmpol package. Hence it is sufficient to include the file 
harmpol.h. Example: 

#- 

#include harmpol.h 
Off statistics; 
.global 

L F = S(R(1,2,1,2) ,N) ; 
#call invmel(S,N,H,x) 
. sort 

On Statistics; 
Print +f; 
. end 

results in 
#- 

Time = 2.72 sec Generated terms = 4 

F Terms in output = 4 

Bytes used = 144 

F = 

- 38/35*z2~3 + l/2*z3~2 - 1/2*H(R(1) ,x)* [1-x] ~-l*z2~2 

- H(R(l,0,l,l,0),x)*[l-x]~-l; 

Of course it should be clear that with these techniques one can do many integrals involving 
combinations of powers of x, denominators with l±x, ln(x), ln(l ± x) and polylogarithms with a 
large variety of arguments. This is usually a soft spot in most computer algebra systems. 

A.4 Color 

This package handles group invariants. The theory behind this is given in ref. j/J. All internal 
variables in this package have a name that starts with the characters cOl in order to avoid potential 
name conflicts. With this package one can, for instance, compute the group invariants that are 
needed in complicated Feynman diagram calculations in such a way that the values for a specific 
representation of a specific group can be substituted at a later stage. This makes calculations more 
general. A simple example of the use of this package is the calculation of a diagram with 14 vertices 
each having three legs, connected in such a way that the smallest loop contains 6 vertices. This 
diagram is unique. If the vertices all belong to the adjoint representation the program would be: 

#include color. h 
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Off Statistics; 
.global 

G gl4 = cOlf (c01il,c01i2,c01i3)*c01f (c01il,c01i4,c01i5) 
*c01f (c01i2,c01i6,c01i7)*c01f (c01i3,c01i8,c01i9) 
*c01f (c01i4,c01il0,c01ill)*c01f (c01i5,c01il2,c01il3) 
*c01f (c01i6,c01il4,c01il5)*c01f (c01i7,c01il6,c01il7) 
*c01f (c01i8,c01il8,c01il9)*c01f (c01i9,c01i20,c01i21) 
*c01f (c01il0,c01i21,c01il5)*c01f (c01il3,c01il9,c01il4) 
*c01f (c01il7,c01ill,c01il8)*c01f (c01il2 , c01il6 , c01i20) ; 

sum cOlil , . . . , c01i21 ; 

. sort 

#call docolor 
. sort 

On Statistics; 
Print +f + S ; 
. end 

and its output is 

Time = 1.93 sec Generated terms = 3 

gl4 Terms in output = 3 

Bytes used = 82 

gl4 = 

+ 1/648*c01NA*c01cA"7 

- 8/15*c01d444 (cOlpAl , c01pA2 , c01pA3) *c01cA 
+ 16/9*cDld644 (cDlpAl , cDlpA2 , cDlpA3) 

This was run on a pentium 300. Note that the customary f has to be given as cOlf etc. Similarly 
the regular Ca is given as cOlcA. The values for the invariants for a given group can be found in 
the paper. As an example we give here the value of E8: 

* 

* Numbers for exceptional algebras 
* 

S [NA+2],eta; 

id c01d644(cQlpAl,c01pA2,c01pA3) = 175/48*c01cA~7*c01NA/ [NA+2] "2 ; 
id c01d444(c01pAl,c01pA2,c01pA3) = c01cA~6*c01NA/ [NA+2] "2 

*(125/27+125/216*c01NA) ; 

* 

* And for E8 
* 

id cOlNA = 248; 
id l/[NA+2] = 1/250; 
id cOlcA = 30*eta; 
Print +f +s; 
. end 

which gives 
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gl4 = 

+ 2075760000*etaT 

with the value of rj usually being two. 
A. 5 Mincer 

This is a pure particle physics package. It concerns the evaluation of massless propagator type 
integrals with up to three loops. In addition there are routines for the expansion in (fixed) Mellin 
moments. Computation time for higher Mellin moments can be rather large. In addition the disk 
requirements increase exponentially with the number of the moment. One needs to specify the 
topology of the diagram and a few other options. A simple example is given by 

#- 

#define SCHEME "0" 
#define T0P0 "la" 
V pl,...,p8,Q; 
#include mincer. h 
Off Statistics; 
.global 

L Dia = 4*p7.p8~2/pl.pl/p2.p2/p3.p3/p4.p4~2/p5.p5/p6.p6/p7.p7/p8.p8*Q.Q; 
. sort 

#call integral ('TOPO') 
. sort 

On Statistics; 
Multiply ep~3; 
id ep = 0; 
print ; 
. end 

which results in 
#- 

Answer in MS-bar 

Time = 0.85 sec Generated terms = 5 

Dia Terms in output = 5 

Bytes used = 64 

Dia = 

- 139/2 - 5/2*ep~-3 - 19/2*ep"-2 - 55/2*ep"-l + 5/2*z3; 

For more details one should have a look at the paper that comes with the package. This mincer 
package has been in use by the author for 9 years now. Over the years it has steadily been improved 
and optimized. The code is basically code for version 2 of FORM. It is this package that was used 
to test the parallel version. 
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